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Transport properties of a superconductor-semiconductor-superconductor (S-Sm-S) junction with 
superlattice structure are investigated. Differential resistance as a function of voltage shows oscil- 
latory behavior under the irradiation of radio-frequency (RF) waves with the specific frequency of 
1.77 GHz regardless of the superconducting materials and the junction geometries. Experimental 
data are quantitatively explained in terms of the coupling of superconducting quasiparticles with 
long-wavelength acoustic phonons indirectly excited by the RF waves. We propose that the strong 
coupling causes the formation of novel composite particles, Andreev polarons. 



Transport properties of Superconductor- 
Scmiconductor-Supcrconductor (S-Sm-S) junctions 
have been studied extensively from both experimental 
and theoretical sides. Even in the S-Sm-S junctions 
where the supercurrent does not flow, diverse phenom- 
ena found in the dynamics of correlated quasiparticles, 
such as multiple Andreev reflections, attract research 
interests [HQ. I n this paper, we report novel composite 
particles, which we call Andreev polarons, formed by the 
strong coupling of quasiparticles with long-wavelength 
phonons that are indirectly excited by radio-frequency 
(RF) waves. Polarons have been recognized so far 
as the composite particles, which are formed in the 
normal conducting (not superconducting) materials or 
dielectric materials by the electron-lattice interaction 
such as the deformation potential due to the acoustic 
phonons and Coulomb interaction in ionic crystals Q • In 
conventional superconducting materials, by contrast, it 
has been unknown that carriers other than Cooper pairs 
can couple with phonons although the phonon-mediated 
attractive interaction between electrons is famous for 
the origin of Cooper pairs. 
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FIG. 1. Schematic cross-sectional view of the sam- 
ples. The superlattice structure consists of two 
GaAs(2 nm)/GaAsNSe(2 nm) layers and three GaAs(l 
nm)/GaAsNSe(l nm) layers. 

Figure [T] shows the schematic cross-sectional view of 



our samples. GaAs/GaAsNSe superlattice (SL) struc- 
tures are sandwiched between an n-GaAs layer and two 
Al superconducting electrodes that form an S-Sm-S junc- 
tion. Each of the SL structure consists of two sets of 
GaAs(2 nm)/GaAsNSe(2 nm) layers and three sets of 
GaAs(l nm)/GaAsNSe(l nm) layers 0]. The carrier den- 
sity and carrier mobility in the n-GaAs layer at room 
temperature are 5xl0 18 cm -3 and 1000 cm 2 /Vs, respec- 
tively. We prepare samples with different lengths of the 
n-GaAs layer (L), which we identify as the length of the 
slit between the two Al electrodes (0.5 /xm and 1.0 /xm), 
while the junction width (W) is speculated to be some- 
what shorter than the designed value of 0.6 /im. We mea- 
sure the differential resistance of the S-Sm-S junction by 
means of standard lock-in technique at the temperature 
of 45 mK under the irradiation of RF waves. Due to the 
impedance mismatch at the handmade loop antenna, the 
irradiated power of RF waves is less than ~1 % of the 
monitored RF sweeper power. The mean free path (£) 
and the coherence length (£n) hi the n-GaAs layer at 
low temperature of 45 mK are estimated to be ~0.1 /im 
and ~0.3 /xm, respectively. 

The measurement results are shown in Fig. [2] Al- 
though we can see some features in the differential re- 
sistance characteristics without RF irradiation, which is 
closely related to the dynamics of superconducting quasi- 
particles, the rates of variability of differential resistance 
in the samples with L =0.5 /urn and 1.0 /xm are 10 % 
and 2.5 % , respectively. These rates of variability are 
significantly small compared with those in ideal S-Sm- 
S junctions, and suggest resistive channel be placed in 
parallel with the junction. Taking into account the fact 
that the mean free path (I) is smaller than the length of 
n-GaAs layer (L), the quasiparticles experience several 
scatterings during the travel in the n-GaAs layer, and 
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FIG. 2. Differential resistances as functions of voltage for 
the length of the n-GaAs layer L —0.5 /im and 1.0 (im. The 
heavy black and thin red lines represent the data without RF 
and with RF irradiation, respectively (Frequency, 1.77 GHz; 
Power, -10 dBm). Voltage is obtained from the measured volt- 
age (V m ) and the measured current (7 m ) using the calibration 
formula V = V m — Rchlm- 



lose the phase information of wave function. Especially, 
in the case where L is longer than ~1 /im, the channel 
resistance of the n-GaAs layer is comparable with the in- 
terfacial resistance at the Al/SL interfaces, and we can- 
not neglect the voltage drop across the n-GaAs layer in 
discussing the quasiparticle dynamics. Therefore, when 
we plot the differential resistance in Fig. [21 we use the 
voltage difference between the two Al/SL interfaces ob- 
tained by the calibration formula; V — V nl — Rchlm- 
Here, V m and I m are the measured voltage and current, 
respectively, and the resistance of the n-GaAs channel 
(i?Ch) is estimated to be SI and 1.4 f2 for the sample 
with L = 0.5 /im and 1.0 /im. Because the resistance of 
the n-GaAs channel (i?ch) is predominantly determined 
by the additional parallel current channel, the procedure 
of our calibration can be validated. 

On the other hand, the coherence length (£/v) is com- 
parable with the length of n-GaAs layer (L), and the cor- 
relation of quasiparticle-pair marginally remains in the 
n-GaAs layer. Therefore, we can discuss the behavior of 
differential resistance characteristics using the Andreev 
bound state formed in the n-GaAs layer. The differ- 
ential resistance characteristics without RF irradiation 
show peaks at V=±V pea k=±110 /aV in the two samples 
shown in Fig. [2] Although this V pea k value is signifi- 
cantly small compared with the twofold superconducting 
gap voltage of bulk aluminum (2Aai/c = 340 /iV), we re- 



gard Vpeak as the twofold superconducting gap voltage of 
Al electrodes (Vpeak = 2A/e), because the quasiparticle 
dynamics drastically change at this value. The energy of 
the Andreev bound state under zero voltage is approx- 
imately equal to the superconducting gap voltage A/e, 
and the peaks/dips appear in the differential resistance 
when the superconducting gap edge of either side of Al 
electrodes coincides with the Andreev bound levels [5j. 

Under the irradiation of RF waves (frequency, 1.77 
GHz; power, -10 dBm), the differential resistance shows 
oscillatory behaviors in both of the two samples as shown 
in Fig. [2] The oscillations are shown only in the voltage 
near Vp ea k= HO /iV, and the voltage period of the oscil- 
lation varies from 25 /iV to 15 /iV, which monotonically 
decreases with increasing \V\. We especially note that 
the minimum voltage period of 15 /iV just coincides with 
twice the energy of the irradiated RF photons (1.77 GHz 
- 7.3 (teV). 

The S-Sm-S junctions without SL structures do not 
show the oscillatory behaviors in differential resistance 
under the irradiation of RF waves, while we obtain al- 
most constant resonant frequency of 1.77 GHz for the 
junctions with SL structures regardless of the supercon- 
ducting materials (Al or Nb [6j) and the junction ge- 
ometries (the length of the n-GaAs layer (L) and the 
junction width (IF)). Therefore, we believe that some 
phonon modes inherent to the SL structure are excited 
by the irradiated RF waves and affect the transport prop- 
erties of the junction. We will discuss later details of our 
scenario that the excited phonons strongly couple with 
the quasiparticles in the S-Sm-S junction and form novel 
composite particles which we call Andreev polarons. 

In Fig. [31 we plot the peak position of the differential 
resistance as functions of RF wave frequency/power for 
the sample with L=0.5 /im. When we slightly change 
the frequency of the irradiated RF waves, the peaks of 
the differential resistance shift together maintaining the 
voltage period almost constant, and take their minima at 
the resonant frequency of ~1.77 GHz. The number of the 
observed peaks is found to decrease as the RF frequency 
deviates from the resonant frequency, and the oscillatory 
behavior totally disappears when the frequency devia- 
tion reaches ^20 MHz. The number of the peaks is also 
changed by RF power. With increasing RF power, new 
peaks appear in both higher and lower voltage regions 
while the voltage period of the peaks is almost constant. 

First we discuss the involved phonon modes. Only 
acoustic modes admit such low-energy phonons that RF 
excitation is possible. Using the sound velocities of acous- 
tic phonons in bulk GaAs (4500 m/s and 2800 m/s for 
longitudinal and transverse mode, respectively, in (001) 
direction Q), the wavelengths are estimated to be 1.6 
/^m and 2.6 /im for longitudinal and transverse mode. 
These wavelengths are comparable with the length of 
the n-GaAs layer (L), and we consider that the standing 
wave of acoustic phonon mode is excited between two Al 
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FIG. 3. Peak position in differential resistances of the sample 
with L=0.5 /im under the irradiation of RF waves, (a) RF 
power = -20 dBm, (b) RF frequency = 1.7694 GHz. 



electrodes by the RF irradiation. However, the acoustic 
phonon modes have poor coupling efficiencies with elec- 
tromagnetic waves because the electric dipole moments 
are not induced in the crystals by the acoustic vibrations. 
Therefore, this excitation needs another phonon mode 
which mediates the RF photon energy to the acoustic 
phonon mode. This mediating phonon mode is required 
to be RF active and inherent to the SL structure. The 
most probable candidate is the interfacial phonon mode 
which runs along the GaAs/GaAsNSe interfaces in the 
SL structure. Due to the elastic boundary conditions 
imposed at the GaAs/GaAsNSe interfaces together with 
the piezoelectricity in GaAs crystals, the interfacial mode 
can be RF active and have the resonant frequency of 1.77 
GHz, which is determined by the sizes of GaAs/GaAsNSe 
interfaces. Although the dispersion relations of the inter- 
facial phonon modes are nonlinear in general, the sound 
velocities in low-energy long- wavelength limit are slightly 
smaller than the transverse sound velocity in bulk GaAs 
layers (i.e. 2800 m/s in (001) direction). Thus, the inter- 
facial phonon mode can couple strongly with the acous- 
tic phonon modes of the bulk n-GaAs layer because the 
phonon energies are close and the pseudo- phase match- 
ing condition can be satisfied over a certain distance. Our 
speculation of the indirect excitation by the mediating 
interfacial phonon modes is also supported by the RF 
frequency dependence of the peak position shown in Fig. 
(3](a). The slight decrease of peak positions in the vicin- 
ity of the resonant frequency suggests the frequency en- 
trainment with a lower-energy phonon mode caused by 
nonlinear phonon-phonon interaction. 



As a result, the superconducting quasiparticles coexist 
with RF-excited long-wavelength acoustic phonons be- 
tween two Al electrodes. Taking the coupling with the 
SL structures into account, both the standing wave of 
the quasiparticle wave function and that of the phonon 
vibration are considered to be limited to the thin sur- 
face region in the n-GaAs layer, and the quasiparticle- 
phonon interaction is greatly enhanced due to the di- 
mensional confinement [8j4l0j. For the formation of the 
novel composite particles, only the coexistence and the 
strong quasiparticle-phonon interaction are essential, and 
the generating mechanism of Andreev polarons is anal- 
ogous to that of standard polarons under the electron- 
phonon interaction [Tll j . which we outlined briefly in the 
following. 

We consider a one-dimensional system of supercon- 
ducting quasiparticles strongly coupled with a specific 
acoustic phonon mode vibration. For the sake of simplic- 
ity, we discuss the transport properties of the junction in 
voltage state using the eigenstates, which exist only un- 
der zero voltage in a rigorous sense. Assuming that quasi- 
particles have discrete energy levels sufficiently separated 
from each other, we apply Born-Oppenheimer adiabatic 
approximation to the energy eigenvalue problem of this 
coupled system, and obtain the following two equations: 
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<p(Q) = E <p(Q). (2) 



Equation ([T]) is Bogoliubov-de Gennes equation with 
quasiparticle-lattice interaction term V{x, Q). For fixed 
value of the phonon coordinate (Q), we can solve Eq. 
([T]) in principle, and obtain the discretized energy eigen- 
value of the quasiparticle system (e(Q)) together with the 
wave functions of quasiparticle eigenstates (u = u(x; Q) 
and v = v(x;Q)). Although we predominantly apply 
the boundary condition of Andreev reflection at the two 
Al/SL interfaces, normal reflection is also admissible due 
to the interface barrier. Equation @ is Schrodinger 
equation for the wave function representing phonon vi- 
bration (<£>(Q)), where harmonic potential is assumed and 
the energy of the quasiparticle system (e(Q)) is included 
in the left hand side. 

We expand the quasiparticle-lattice interaction term 
V(x,Q), together with e(Q), u(x;Q) and v(x;Q), into 
the power series of normalized phonon coordinate Q/A 
where A = a(hfMQ) 1 ^ 2 is the characteristic length which 
corresponds to the amplitude of the zero-point vibration 
in the phonon mode with numerical multiple a > 10; i.e. 
V(x, Q) ee Vx(x)(Q/A) + • • • , e(Q) = e + e x (Q/A) + 
(u,v) ee (u ,v ) + («i,«i)(Q/A) + • • •, respec- 
tively. It is noteworthy that the zeroth-order term of the 
quasiparticle-lattice interaction is omitted (Vq(x) ee 0) 
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because it has been already included in H and A in 
Eq. (pj. Using the standard perturbation technique, 
we obtain the energy of the quasiparticle system as 
e(<9) ~ eo - (Vi)Q/A + 0{Q 2 /A 2 ). Here we note that 
V(x, Q) is a 2x2 matrix in general, and define the expec- 
tation value (Vi) as (Vi) = J dx (uq, Vq)Vi(x)(uo, vo) t - 
We can diagonalize the Hamiltonian in Eq. ([2]) with 
e(Q) « e - (Vi)Q/A by putting Q> = Q - (Vi)/MCl 2 A, 
and obtain the energy eigenvalue of the total system i.e. 
the energy of Andreev polaron as 

^ = e*-^ + ttl(n+±). (3) 

Here the quantum number n can be considered as the 
number of the phonons which one quasiparticle is dressed 
with. 
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FIG. 4. Schematic diagram of Andreev polarons together with 
the adiabatic potential (parabolic line). The superconducting 
gaps (A) in the Al electrodes are also shown. 



Thus, with each voltage period of 5V — 2Ml/e, the dif- 
ferential resistance characteristics of the junction is mod- 
ulated because of the novel resonant condition. We can 
state that the new channel concerned with the n-th An- 
dreev polaron level participates to the transport of the 
junction with each voltage period of SV. 

Because the stabilization energy of Andreev polaron is 
larger than the bare energy of the quasiparticle eo (~ A) 
and the zero-point energy of the excited acoustic phonon 
Ml/2, the first term in the right hand side of Eq. (jj) is 
negative. From the peak position with the largest |V|, 
which corresponds to 2Eo/e, we estimate the magnitude 
of quasiparticle-lattice interaction (Vi)(Q)/A as large as 
~0.4 meV. When the irradiated RF power is increased, 
the expectation value of the phonon coordinate (Q) is 
increased, which is consistent with the appearance of new 
peaks shown in higher and lower voltage regions (Fig. 

EDO). 

In conclusion, we investigate transport properties of a 
superconductor-semiconductor-superconductor (S-Sm-S) 
junction with superlattice structure. The differential re- 
sistance is measured as a function of bias voltage which is 
found to oscillate under the irradiation of radio- frequency 
(RF) waves with a specific frequency of 1.77 GHz. This 
oscillatory behavior is quantitatively explained by the 
strong coupling of superconducting quasiparticles with 
long-wavelength acoustic phonon mode vibration which 
is indirectly excited by RF waves. We propose the re- 
sultant formation of novel composite particles, Andreev 
polarons. 

We thank Dr. D. Tsuya and Dr. E. Watanabe in Na- 
tional Institute for Materials Science (NIMS) Nanotech- 
nology Innovation Center for nanofabrication. 



In Fig. SI we show the schematic energy diagram of 
Andreev polarons under zero voltage together with the 
adiabatic potential (Mil 2 Q 2 /2 + e(Q)) which the phonon 
mode vibration is subjected to. Due to the quasiparticle- 
lattice interaction, the adiabatic potential takes its mini- 
mum at a finite phonon coordinate of (Q) = (Vi) /MVL 2 A, 
which is lower than the bare energy of quasiparticles 
(eo) by (Vi) 2 /2M7. The energy-level spacing of Andreev 
polarons is approximately equal to that of the involved 
phonon ft£l to an accuracy of 0(Q/A). 

In order to describe the transport properties of the 
junction under the voltage state, we have to obtain the 
occupation probabilities of the involved Andreev bound 
states, which are different from the equilibrium Fermi 
distributions in general, and calculated by Boltzmann 
equation self-consistently, for example [3, 0|. However, 
the resonant voltages where peaks/dips are shown in dif- 
ferential resistance characteristics can be obtained as 



eV = 2E„ 




(4) 
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